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ABSTRACT 

Euler's identity is shown to give a relation between the zeros of 
the Riemann-zeta function and the prime numbers in terms of $(a;) = 



_i 

X 4 



2^i£e- p2 ™ln(j») - 1 



on x > 0, where the sum is over all primes p. If $ 
is bounded on x > 0, then the Riemann hypothesis is true or there are infinitely 
many zeros Re Zk > |. The first 21 zeros give rise to asymptotic harmonic 
behavior in $ defined by the prime numbers up to one trillion. 



1. Introduction 

The Riemann-zeta function £(z) is a complex function of z defined in Re(z) > 1 by the 

C(z) = 1 + — + — + •■■ (1) 

M ' 2 3 

over the integers. The zeros Zk of ((z) comprise the negative even integers, —2, —4, • • •, and 
an infinite distribution of nontrivial zeros z^ = + iyu in the strip < < 1. 



sum 



A general approach to find zeros of as endpoints of continuation (iKellerl 119871 ). Here, 
we set out to analyze the Zk of as endpoints of continuation along a trajectory z(X), whose 
tangent z'(X) satisfies 

in response to a choice of initial data z(0) = zq, is a zero of £(z). The trajectory z(X) 
continues whenever r ^ 0, i.e., whenever Q(z) = exists and is finite. Continuation 

halts only at a zero of l/0(z), corresponding to a zero of C( z )i au °f which are isolated. 
Continuation (j2J) can equivalently be formulated in terms of the prime numbers by 



£0) = E\n(p)p- z (Re z > 1), 



(3) 
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representing the leading term the analytic extension of the sum 

mmz) = = ~ C M (4) 

by virtue Euler's identity 

C(z) = n (1 + P - z + P - 2z + •■■)= n (i - v - z y x . (5) 

All known nontrivial zeros determined by numerical root finding satisfy Re z\~ = | to 
within numerical precision, the first three of which are z\ — | ± 14.1347z, z 2 = | ± 21.0220z, 
z 3 = | ± 25.0109i. By the symmetry 

C ( s ) = C(l-s)^p xW=T»r(|), T{z) = t*~ l e-Ht, (6) 

it suffices to study zeros in the half plane Re(z) > \. Fig. CD illustrates root finding by (T5]) 
for the first few zeros. 

In this paper, we study £(z) in ([3]) in the relations between the distribution of primes 
and the location of Zk by Euler's identity in the expanded form 

£(z) = -^4-£ m>2 £(mz). (7) 

The results will be formulated in terms of a density and normalized density function of the 
primes defined on the positive real numbers x > 0, given by 

$(x) = 2v/i0( f ) ' 1 , 0(x) = Se-^logb), (8) 

where summation in <p{x) is over all primes. The function has the property that it 
possesses poles at the zeros z k of ((z). In particular, the pole at z = 1 is a c onsequence 
of the Prime Number Theorem applied to the Chebyshev function JPusartlll999l ). by which 



2^Jx(f)(x) — 1 = o(l) as x approaches zero. In our application of a Fourier transform to 
the Zk introduce a sum of residues 

Z{\) = S Tfc e- A ^-5) (9) 

defined for all A < 0, where the 

ik = i{zk), -y(.z) = -¥ l , (io) 

7T2 

are absolutely summable by Stirling's formula and the asymptotic distribution of z^. 
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Fig. 1. — Shown are the trajectories of continuation z(X) in the complex plane z by numerical 
integration of §2§ with initial data z = 1 + ni (n = 1, 2, 3, • • •) indicated by small dots on 
Re(z)=l. Continuation produces roots indicated by open circles, defined by finite endpoints 
of z(X) in the limit as A approaches infinity. The roots produced by the choice of initial data 
are the first three on Re z = ~ and -2 and -4 of the trivial roots. 



Theorem 1.1. In the limit as x > becomes small, we have the asymptotic behavior 



2vri 



Z (In VS) 



+ 



37T6 



X 12 



O \ X 12 



TT 



Corollary 1.2. If $ is bounded, then the Riemann hypothesis is true or there are 
infinitely many zeros Re Zk > h. 

A similar relation bet ween the distribution of Z\. and the primes is given by (IHadamard 
18931 : Ivon Mangoldtl fl895h 



ip(u) — u ^u Zk 2 ln(27r) + In y/l — u~ 2 



'U Z k y/U 

based on the Chebyshev functions 

i)(u) = S p fc< u ln(p), = S p < u lnp, 



121 



(13) 
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where the sum is over all primes p and integers k. In Theorem 1.1, normalization of $(x) in 
Theorem 1.1 is by x * according to (jSJ) and Z is absolutely convergent for all x > 0, whereas 



in ( fl2l) normalization is by \/u and the sum £2-=-^ is not absolutely convergent. Similar to 



Corollary 1.2, it will be appreciated that the left hand side of f fl2|) will be bounded in the 
limit of large u if the Riemann hypothesis is true. 

In §2, we give some background on (,{z). in §3 we introduce an integral representation of 
£(z) and derive one of its integral properties associated with the singularity of ((z) at z — 1. 
In §4, we apply Cauchy's integral formula to ((z) to derive a sum of residues associated by 
the Zk- The proof Theorem 1.1 is given in §5 by a Fourier transform and asymptotic analysis 
of the expanded Euler's identity. In §6, we report on a direct evaluation of $(x) using the 
primes up to one trillion, to show the harmonic behavior as as encoded in Z by the first few 
zeros Zf.. We summarize our findings in §7. 



in 



2. Background 

Our analysis use s some w e ll known properties o f C(z) and its z e ros Zu = g±+ iyk, e.g. 



Titchmarshl (I1986T ): iLehmeri (119881 ): iDusartl (119991 ): Keiperl fll992f ); iFordl fl2002f )l 



Riemann provided an analytic extension of £ (z) to the entire complex plane by expressing 
each term n~ z in terms of T (|) , 



CO*) 



7T 2 



r(!) 



where 



x 2 1 6i(x)dx, 



(14) 



9 1 (x) = 9 -^-^, 6{x) 



(15) 



satisfies the asymptotic property 9i(x) ~ as x approaches zero by the identity Oi^x^ 1 ) = 
tJx6(x) for the Jacobi function #(x)o Evaluation of the integral (I14p on Re z > 1 hereby 
reveals the meromorphic expression (e.g. Borwein et al. 2006) 



C(z) 



) = i% \w^) +f{z) )' f{z) = l i x "~ 1 + x ~^) e ^ dx 



(16) 



When z — n is an integer, the pre- factor 3 - % will be recognized to be half the surface area S(n) of a 
unit n— sphere. 
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representing a maximal analytic continuation of ((z), and the presence of a simple pole at 
z — 1 with residue 1. 

Riemann also introduced the symmetric form Q(z)((z), Q(z) = \z{z — l)7r~"ir (|) 
satisfying Q(z)((z) = Q{1 — z)((l — z ), whereby 



((z) = n*-\(l-z)- 12 2 ' 



nX(l-z) 



r( 



n*-i2 z -\(l - z) 
cosQttz) r(f)r(i + f) ~ cos (|ttz) f£z) 



(17) 



using T(| - |)r(| + §) = tt^-t and r(z)r(z + \) = 2 l - 2 \/^T (2z). Along the vertical line 
z = 1-iy, C(z) ^ (|Littlewoodl1l922l . Il924l . Il927l : Iwintner!ll94ll ) . which allows consideration 
of its logarithmic derivative, 

^4 = - C ?~ Z ] - + - tan ( _ ^( Z ) (18) 

C(*) C(l -2) 2 V 2 / ^ v ; v ; 



in terms of the digamma function 



T\z) 
T(z) 



ln(z) + 0{z~ l ) 



(19) 



in the limit of large \z\. 

Lemma 2.1. In the limit of large y, the logarithmic derivative of ((z) satisfies 

C'iiy) C'(l - w) 



City) C(l - »y) 



0(lBJ/). 



(20) 



Proof. The result follows from (Tl9|) and (I18p . □ 

Lemma 2.2. Along the line z = iy, we have the asymptotic expansion \^f{iy)\ 
-a"^ y in the limit of large y, whereby the 7^ are absolutely summable. 



Proof. Recall fflUj) and the asymptotic expansion T(z) = y/2nz z ze z [1 + 0(z r )] with 
a branch cut along the negative real axis. In the limit of large yk, and the asymptotic behavior 



of the zeros Zk, whereby yk ~ f 21 !, and hence 



e Ink 



where we used | arg Zk \ 



~ f in the 



limit of large k. It follows that the 7^ are absolutely summable. Numerically, their sum is 
rather small, £|7fc| = 3.5 x 10~ 5 based on a large number of known zeros z^. □ 



Lemma 2.3. In the limit of large y, we have 



i{w) 



C(iy) 



O (y *e 2^ lnyj . 



(21) 
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Proof. By Lemma 2.1-2, we have 



liiy) 



C'(iy) 



C(iy) 



< h(iy)\ 



C'(i - iy) 



C(i-iy) 



C'(i-«v) 



C(i-*y) 



+ 0(lny) (22) 



asymptotically for large y. By ( jRichertlll9671 ; IChengill999l ; lTitchmarshlll986l ; lBroughanll2005l ) 

CO—iy) 



C0--iy) 

on y > 5 for some positive constants c, 5. □ 

These observations give rise to the following. 
Proposition 2.4. The Fourier transform 



< c(ln7/) 3 (lnlny)' 



(23) 



L(A) = PV 



7(%)g%"^ + ln(2.) 



(24) 



where PV refers to the Principle Value about y = 0, is well defined and F(X) — > in the 
limit of large A. 



Proof. The result follows from Lemma 2.1-3 and the Riemann-Lebesque lemma. □ 



3. An integral representation of 

We begin with an integral representation of in terms of (p(x) and given in 
OH]), following steps very similar to those by Riemann leading to ( II 6p . To this end, we define 
g(z) = g 1 (z)+g 2 (z), 



/U /-oo 
$(e 2A )e A( ^dA, = 2 / 0(e 2A )e Az dA (25) 

■oo JO 



for Re z > 1 and, respectively, all z. 

Proposition 3.1. On Re 2; > 1, we have 

S ln(p)p- 2 = JLL. I _L_ + g ( z ) I . (26) 

2 

Proof. Following steps similar to those leading to f fT6|) . we have £(z) = -^y / °° x%~ l <fi(x)dx, 
that may be decomposed in integration over (0, 1] and [1, 00), 

9i(z) = \ I x^®{x)—, g 2 {z) = [ x*~ x <i>{x)dx, (27) 



2 In y ' x 
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where gi{z) includes regularization about z = 1 made explicit in fl26|) by -j^-, where the 
singularity at z = 1 in £ can be seen to result from the Prime Number Theore m in the form 



of 2yfx(j)(x) — 1 = o(l) on the basis of the asymptotic behavior of the function (jDusartlll999l ). 
The substitution x = e 2A gives ( 1251) . □. 

Lemma 3.2. TTie analytic extension of gi(z) extends to the real values z > ~. 

Proof. In a neighborhood of < z < 1, we may write 

C(z) vrf 



C(z) r(|) 



i 



+ uAz) 



(28) 



where 1*1(2;) is analytic about z > 0. With 2 = a + ib, the second term on the right hand 
side in the expanded Euler's identity ([7]) satisfies 

S m > 3 l£(m^)l < E n > 3 1 _ n _ a < — - , (29) 

whereby it is bounded in Re z = a > \. Since the second term C(2#) in (E]) is analytic in 
Re z = a > |, it follows that g(a) as defined in Proposition 3.1 is analytic on a > h. In view 
of the analytic and finite behavior of the right hand side in (j7]), the first and second term on 
the left hand side in ([7]) remain balanced as a approaches | from the right, giving 

£(a) = -— L- + U2 (a), (30) 
2a — 1 

where «2(a) is analytic at a = \. As a approaches | from the right, we have 

0i (a) = 1 : +u 3 {a), (31) 
2a — 1 

where 7/3(0) is analytic about a = |. □ 

By Proposition 3.2, the behavior of the sign of $(x) in a neighborhood of x = is 
relevant. If $(x) is of one sign in some neighborhood, then gi(z) has an analytic extension 
into Re z > \ with no singularities, implying the absence of z& here. However, this requires 
information on the point wise behavior of $(x), which goes beyond the relatively weaker 
integrability property ( 13T|) . To make a step in this direction, we next apply a linear transform 
to (JTj) to derive the asymptotic behavior of $(x) in terms of the distribution Zk- 



4. A sum of residues Z associated with the non-trivial zeros 

We next define 

M2) .£li)ai) + _L (32) 

7T2 Q{Z) Z-l 
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and its Fourier transform 



H(\) 



a+ioc j 



(33) 



Lemma 4.1. The function h(z) has a simple pole at z = 1 with residue 1 and simple 
poles at each of the nontrivial zeros of ((z) with residue 7^. 



Proof. As a consequence of (e.g. Borwein et al. 2006) 



Uz)= l -z{z-l) V -^-Q{z\ ff\ 

2 7T2 Cl(^) 



1 1 
+ - I • 

z- z k z k 



(34) 



where B is a constant, so that 



r (!) CO*) + _2. 



7T2 



Here 



C(*) 



A 



r(I) 



7T2 



1 1 
+ — 

z — z k z k 



+ A. 



(35) 



1 L r(f)\ r(|) 2 r(|) 

2 — 1 \ 7T2 / Z7ri 7rfln7r 



(36) 



includes contributions from the logarithmic derivative of the factor to £(2) in (134p . whose 
singularities are restricted to the trivial zeros of ((z). □ 

Proposition 4.2. The Fourier transform of h(z) over Re z > supa^ satisfies 

H(\) = e-?Z(\) + 0(1) (37) 

in the limit of large A < 0. 

Proof. Consider contour integration over the rectangle Re z = 0, Re z = a and z = x±iY 
(0 < x < a). The line integral over z = satisfies O ( e ^) by Proposition 2.4. Integration 
over z = x + iY (0 < x < a) gives 



riY+a j 

e~* Y / h{z)e- Xx — 
JiY 27n 



e - lY r Y+a mm^ +0 (Y-i) 

U vrf ((z)2m +U[Y > 



(38) 



where we choose Y to be between two consecutive values of y k . We have 



— r (.)«,)* . 7i f + «W^ ln(2a _ i; 



f*y ^ C(z) 2 ™ 



iY C(z) 2m 2vn 



4a% fc - Y) ; 

1 H h TTi 

I -2a 



(39) 
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In the limit as k approaches infinity, y k — Y approaches zero and \j k \ becomes small by 
Lemma 2.2., whence 

( f Y+a ~ r ¥+a ) T M^\^. ~ H2a - l)Im lk = O ( ln(2a - l)^e^A . (40) 

\JiY J-iY J 7T 2 C{Z) 2m V V Vk J 

The result now follows in the limit as k approaches infinity, taking into account the residue 
sum e"2Z(A) associated with the z k and absolute summability of the 7^. □ 



5. Proof of Theorem 1.1 

Following the expanded Euler identity ©, we define 

r ( z -) 

f(z) = -g 2 (z) - h(z) - -^Lv) m > 2 i(mz) (41) 

7T2 

for all z away from the z k , its Fourier transform 

ra+ioo . 

F(X) = / Me^— (42) 
on a > supafc, and the coefficients 

c m ( z ) = ^4 , C m = m-^-^T ( — ) . (43) 

r(i) 

for all m > 1. We remark that in particular, C 2 = — H^- and m~ 1 c m (m~ 1 z) = 1 + 

(I In 7r + I7) z + O (z 2 ) has a well defined limit and C m — > 2 in the limit as m becomes 
arbitrarily large. 

Lemma 5.1. TTie decomposition of the sum £ m >i£(mz) = £(z) + S m >2^(m2) is we//- 
defined on Re z > \. 

Proof. By the Prime Number Theorem, p k ~ fcln(fc), whereby summation over the tails 
k > n defined by large n satisfy 



ln(p^ 



J k>n 2a k>n 
Pk 



1 lnln(&) 



fc 2a ln(fc) 2a - 1 (k\n(k)) 



2d 



< 00 (44) 



whenever a > |. Hence, for 2 = a + iy, |Sln(p)p 2 ^| < £ ln(p)p 2a < 00 whenever a > 
f. It follows that |S m > 2 e(m^)| < S m > 2 S p p-( m - 2 ) a ln(p)p- 2a < £ m > 2~ m £ p ln(p)p- 2a = 
T> p \n(p)p~ 2a < 00 on Re z > ^, whereby the sum £ m > 2 £(mz) is well-defined. Given the 
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analytic continuation of S m >i^(m2) to Re z > | by the logarithmic derivative of the (analytic 
continuation of) ({z), it follows that the (analytic continuation of) £(z) is well-defined. □ 

Lemma 5.2. F(X) = e"^$(e 2A ). 

Proof. By Proposition 3.1, the Fourier transform 

a+ioo j 

9i{^— (45) 

-oo 2m 

exists uniquely on Re z > 1. Here it equals F(X), since f(z) is the analytic extension in 
Re z > | of in Re z > 1. □ 

Lemma 5.3. For any m > 2, £/ie Fourier transform of Cto ^ over /?e 2; = a > ~ satisfies 

DM) = C m e-^ + o(l) (46) 

Proof. We have 

A.(A) = «* f J^e-"^ + CU^tt-i). (47) 

The integral ()47l) exists by virtue of a removable singularity of c m (z) at 2 = and decays 
along the imaginary axis when m > 2 by the Riemann-Lebesque lemma. □ 



We now consider the Fourier transform of the expanded Euler identity in the form (|41 
over Re z = a on a term by term basis. To organize terms, we define 



whereby 



r (2) 

r N (z) = g 2 (z) + Y, m > 2 c m (z)g(mz) H ^-T, m > N+1 ^(mz), (4£ 

7T2 



/(*) + sILi-^zt + = ( 49 ) 



and the Fourier transform 

/•a+ioo j 

^tv(A) = / r N {z)e' Xz — . (50) 
Lemma 5.4. R^[X) = o(l) m the limit of large A < 0. 

Proof. Since r^{z) is analytic in Re z > |, its transform i?jv(A) exists and is uniquely 
defined for a > \. The result follows from the Riemann-Lebesque lemma. □ 
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Proof of Theorem 1.1. The Fourier transform of the expanded Euler identity in the form 
(S9D is 

F(X) + D 2 {X) + D 3 (X) + R 3 (X) + H(X) = 0. (51) 
By Propositions 5.1, 4.2 and Lemmas 5.1-5.4, we have 

e -|$( e 2A ) + C 2 e-^ + C 3 e-$ + e~^Z(X) = 0(1) (52) 
With x = e 2A , Theorem 1.1 now follows. □ 

6. Numerical illustration using the primes up to one trillion 

Fig. [2] shows the first harmonic in the normalized density <&(x) for small x evaluated 
for the 37.6 billion primes up to y — 10 12 , allowing x down to 2.6 x 10 -23 (A = —26) in 
view of the requirement for an accurate truncation in 4>(x) as defined by jSJ). The harmonic 
behavior is emergent in 

R(x) = $(x) -C 2 -C 3 x^. (53) 

To search for higher harmonics Z^X) associated with the zeros Z{ in Ae[— 26, —11.7759] = 
[Ai — A2, Ai + A2], we compare the spectrum of $ (e 2A ) by taking a Fast Fourier Transform 
with respect to a, 

X(a) = Ai + A 2 cosa (ae[0, 2ir]), (54) 

and compare the results with an analytic expression for the Fourier coefficients of the Zi(X) 
(i = 1,2,..-), 

c ni [X u A 2 ] = 2Re{(-t) n ll e- a ^J n (-X 2 z t )}, (55) 

where J n (z) denotes the Bessel function of the first of order n. Fig. [3] shows the first 21 
harmonics in our evaluation of $(x), which appears to be about the maximum that can be 
calculated by direct summation in quad precision. 

7. Conclusions 

The zeros Zk = ak + iyk of the Riemann-zeta function are endpoints of continuation: 
i(zk) = iff \/Q(z k ) = 0. By Euler's identity ([5]), continuation and hence the endpoints are 
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controlled by an integral of a regularized sum $(e 2A ) over the prime numbers as defined 
by ®. 

The contributions of the zeros Zk introduce asymptotic harmonic behavior in $ ( e2A ) 
as a function of A < as defined by the sum Z(\) of residues associated with the z k , 
shown in Figs. 1-2. This harmonic behavior in the prime number distribution appears to be 
exceptionally dense: primes up to 4 billion are needed to identify the first 4 harmonics, up 
to 70 billion for the 10 and up to 1 trillion for the first 21. It appears that the prime number 
range scales exponentially with the number of harmonics it contains. Conversely, knowledge 
of a relatively small number of zeros z k provides information on exponentially large primes 
numbers. 

Theorem 1.1 introduces a correlation between $, defined over the primes, and Z, defined 
over the z^. Suppose there are a finite number of zeros z k in Re z > \. We may then consider 
k* for which a k * = maxa^ gives rise to dominant exponential growth in Z(X) in the limit 
as A < becomes large. This observation leads to Corollary 1.2. Z can remain bounded in 
x > only if the Riemann hypothesis is true, or if Z(X) remains fortuitously bounded as an 
infinite sum over a k > \ with no maximum in a < 1. 
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Fig. 2. — The top window shows the normalized density $ (e 2A ) on Ae[— 26, —11.7756] and 
its leading order approximation 1.3616 + 1.5332ee. The asymptotic harmonic behavior is 
apparent in the residual difference (153]) between the two, shown in the bottom two windows, 
including the period of 2.2496 in A associated with the first zero z* — ~ ± 14.1347z. 
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index n 

Fig. 3.— Shown are the absolute values of the Fourier coefficients c n [Ai,A2] of $ (e 2A ) 
obtained by a Fast Fourier Transform (FFT) of (1531) on the computational domain ( l54lh 
where Ai = —26, A 2 = —11.7756] covers 32 periods of -Zi(A) (dots), on the basis of the 
37,607,912,2019 primes up to 1,000,000,000,0039. The resulting spectrum is compared 
with the exact spectra c n j[Ai,A 2 ] of the Zi(\) given by the analytic expression ( 1551) for 
i — 1, 2, 3, • • • (continuous line). Shown are also the individual spectra of Zi(X) for i — 1,8 
and 15 associated with the zeros z\, z§ an d z 15 . The match between the computed and 
exact spectra accurately identifies the first 21 harmonics of Z(X) in $ out of 22 shown, 
corresponding to the first 21 nontrivial zeros Zi of ((z). 



